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Abstract 


We report on the results of three observations of the binary X-ray pulsar SMC X-l with the 
Ginga satellite. Timing analyses of the 0.7 1-s X-ray pulsations yield Doppler delay curves 
which, in turn, provide the most accurate determination of the SMC X-l orbital parameters 
available to date. The orbital phase of the 3.9-day orbit is determined in 1987 May, 1988 
August, and 1989 August with accuracies of 11 s, 1 s, and 3.5 s, respectively. These phases 
are combined with two previous determinations of the orbital phase to yield the rate of change 
in the orbital period: P orb /P orb = (-3-34 ± 0.023) x 10' 6 yr 1 . An interpretation of this 
measurement and the known decay rate for the orbit of Cen X-3 is made in the context of tidal 
evolution. Finally, a discussion is presented of the relation among the stellar evolution, orbital 
decay, and neutron-star spinup time scales for the SMC X-l system. 
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1. Introduction 

Our observations of SMC X-l were motivated by the similarity of this binary X-ray 
system with the Cen X-3 system, which led us to speculate that the orbital period of the former 
system may be changing relatively rapidly as is the case for Cen X-3 (Porb^orb = 'i-8 x 10 ^ 
yr l : Kelley et al. 1983b). The time scale for orbital decay of 5 x 10 5 yr in Cen X-3 is quite 
short by astronomical standards. Studies of such rapid changes in orbital period are important 
to understand (i) the underlying cause, e.g., tidal torques related to the Darwin effect (Darwin 
1879), and (ii) the lifetime of massive X-ray binaries as the Roche lobe moves through the 
stellar atmosphere (only -2000 years in the case of Cen X-3). In particular, we are interested 
to determine if such rapid changes in orbital period are universal among massive X-ray 
binaries. 

SMC X-l is one of only three detected extragalactic binary X-ray pulsars (along with 
LMC X-4 and A0538-66; Primini, Rappaport, and Joss 1977, Bonnet-Bidaud and van der Klis 
1981, and references therein). The optical counterpart has been identified with a BO I 
supergiant (Sk 160; Webster et al. 1972, Liller 1973). X-ray eclipses which occur every 3.89 
days were discovered with the Uhuru satellite by Schreier et al. (1972). SMC X-l exhibits 
low-intensity states in addition to high-intensity states (e.g., Seward and Mitchell 1981, 
Bonnet-Bidaud and van der Klis 1981); Gruber and Rothschild (1984) suggested, on the basis 
of HEAO 1 (A4) data, that the source has a 60 day long-term periodicity similar those found in 
LMC X-4 and Her X-l. 

The source brightness, together with the short pulse period (0.71s; Lucke et al. 1976) 
have made SMC X-l an excellent object for pulse-timing studies. Such studies were initiated 
with SAS-3 and combined with measurements of the optical velocity curve to yield a projected 
orbital radius a x sin i = 53.46 i 0.05 ^t-s, neutron star mass M x = 0.8 - 1.8 Mg, companion 
mass M c - 19 Mg, and companion radius Re - 18 Rg (Primini, Rappaport, and Joss 1977). 
A list of orbital phase zero measurements for SMC X-l has been given by Bonnet-Bidaud and 
van der Klis (1981); this list contains results from observations with Uhuru (in 1971), 
Copernicus (1975), COS-B (1976), and Ariel 5 (1976), in addition to the SAS-3 measurement. 
The first three of these orbital phase determinations were based on eclipse data, while the latter 
two were derived from more accurate pulse-timing analyses. These orbital phase 
determinations are listed in Table 1. Subsequently, additional observations of SMC X-l were 
conducted with HEAO 1 (Gruber and Rothschild 1984), the Einstein Observatory (Darbro et 
al. 1981), and EXOSAT (the EXOSAT log). To our knowledge, however, no additional 
orbital phase determinations have been reported from these latter observations. 
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Our Ginga observations of SMC X-l are described in §11 of this paper. Analyses and 
results pertaining to the timing data and pulse profiles are presented in §111. The results are 
discussed in §IY. 

2. Observations 

The observations of SMC X- 1 were conducted with the Large Area Counters (LAC) on 
the Ginga satellite (Makino et al. 1987, Turner et al. 1989) during three epochs: 3 days in 1987 
May, 8.4 days in 1988 August-September, and 4 days during 1989 July-August. Due to earth 
occultations, crossings of the South Atlantic Anomaly and other high background regions, and 
performance of telemetry, attitude control, and sky monitoring functions, data on SMC X-l 
were obtained with roughly a 30% duty cycle. 

The 3-day 1987 May observation was performed during the initial operational test phase of 
the Ginga satellite, and the data useful for pulse timing were obtained during contact orbits 
(i.e., the ~5 consecutive orbits per day in which the satellite passes over the Kagoshima, Japan 
ground station). Much of these data were taken in MPC-3 mode (12 energy channels with 7.8 
ms time resolution), but the data obtained during the -10 minute intervals of actual 
satellite-ground communication were usually taken in MPC-2 mode (48 energy channels with 
62.5 ms time resolution). For the pulse-timing analysis reported here the MPC-3 data were 
converted to 62.5 ms time resolution by combining time bins, and the energy channels of both 
modes were combined to yield count rates in the 1-7, 7-14, and 14-23 keV bands. During this 
observation the X-ray intensity of SMC X-l was approximately constant at -90 counts s' 1 
(2-37 keV; - 8 joJy at 5 keV). 

The 1988 August-September observation was specifically tailored to monitor the X-ray 
flux of SMC X-l as continuously as possible, while obtaining data useful for pulse timing. 
Data were obtained during the intervals 1988 August 24.9 - 28.3 (UT) and 1988 August 31.1 - 
September 2.7 (UT), but the second interval included an X-ray eclipse. During remote orbits 
(i.e., the —10 consecutive orbits each day in which the satellite does not pass over the 
Kagoshima ground station) the data were obtained exclusively in PC mode (2 energy 
channels). The energy limits for these channels were chosen to be 1-6 keV (with 31.2 ms time 
resolution) and 6-24 keV (with 62.5 ms time resolution). [John: which channels do these 
correspond to?] During contact orbits the observations were similar to those in 1987 May, 
except that MPC-3 data were generally taken with 62.5 ms time resolution. The X-ray intensity 
of SMC X-l outside of eclipse varied between -25 and -300 counts s' 1 (2-37 keV; 2 - 25 |iJy 
at 5 keV). The intensity was particularly low during the time intervals MJD 47398.4 to 
47398.7 and 47400.5 to 47401.3 (MJD = JD - 2,400,000.5). 

Most of the 1989 data were taken with the MPC 2H and MPC 3M modes which provide 
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counts in 62.5 ms time bins. Additional data were accumulated in the MPC 2M, MPC 3L, or 
MPC 3H mode which provided time resolutions of 8 to 500 ms and either 12 or 48 pulse 
height channels. Only those data having a time resolution of 7 .8 or 62.5 ms have been utilized 
for our pulse-timing analyses. Furthermore, these data were rebinned into 10 energy channels. 
The equivalent X-ray energies corresponding to these 10 pulse height channels are given in 
Table 2. The X-ray intensity of SMC X-l was approximately constant at -300 counts s 1 
(2-37 keV; - 25 pJy at 5 keV) for most of the observation outside of the eclipse. 

3. Analysis and Results 

3.1 Temporal Analysis 

Pulse timing analyses have been performed with non-background subtracted data which 
have been selected to exclude the times of the X-ray eclipses, all instances in which the 
detectors were not operating stably, and all readily apparent noise spikes. 

The pulse-timing analysis basically proceeded as follows. The time of observation of each 
bin of data was converted to coordinated universal time (UTC) and corrected for the relative 
motion of Ginga and the pulsar according to the known satellite orbit and a provisional orbital 
ephemeris for SMC X-l. This yields the times of observations that would have been recorded 
by an instrument moving with an approximately constant velocity relative to the pulsar. The 
resulting observation times were then used to fold subsets of the data according to a provisional 
pulsar rotation period and period derivative. Different pulsar orbital and rotational ephemerides 
were used to reduce the data from each of the three observations. Details of this reduction 
procedure follow. 

The observation times were converted from satellite on-board clock time to UTC. The 
relation between satellite on-board clock time and UTC was established to an accuracy of 1 ms 
by comparison of the satellite time with a ground-based standard clock once during each of the 
five times that the satellite passed over the Kagoshima ground station each day. The stability of 
the Ginga on-board clock (-1 x 10‘ 6 ) was sufficient to maintain 1 ms accuracy over the 
100-minute interval between consecutive ground contacts (Deeter and Inoue 1990). For remote 
orbits, however, the use of a mean clock rate between successive clock comparisons (separated 
by some 17 hours) can result in errors of up to -20 ms in the conversion of satellite time to 
UTC. Variations in the Ginga clock rate are largely due to variations in the internal temperature 
of the satellite (Deeter and Inoue 1990). Using the temperature housekeeping data, we were 
able to estimate the clock rate and thereby improve the accuracy of time assignments to -2 ms. 

The observation times were corrected for satellite motion relative to the center of the earth 
using the Ginga orbital ephemeris, and then transformed to the solar system barycenter using 
the CfA planetary ephemeris (Chandler 1991) or the JPL planetary ephemeris (reference). In 
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making these corrections we used the coordinates of Sk 160 reported by Clark et al. (1978). 

The observation times were also corrected for a provisional pulsar orbit. The results were 
then used to fold subsets of data from time intervals typically -2000 s in duration according to 
a provisional pulse ephemeris. The provisional parameters were sufficiently precise to fold a 
—2000 s interval of data with smearing of < 0.1 pulse periods. For the 1987, 88, & 89 
observations, a total of 18, 149, and 61 pulse profiles were constructed, respectively. 

For each of the three observations a pulse template was constructed by averaging the 
individual profiles (1987 and 1988 data) or simply by utilizing one of the pulse profiles having 
particularly good counting statistics (1989 data). Each template was then cross correlated with 
the individual profiles from that observation to obtain pulse arrival-time delays. The resultant 
Doppler delay data (with the delays due to the best-fit orbit added back in) are plotted in Figure 
1 . 

We fit the pulse arrival time delays to circular and eccentric orbit models to obtain 
corrections to the provisional orbital and pulse period ephemerides. All of the models included 
a separate parameter for the mid-eclipse (T^) epoch for each of the 3 observations as well as 
parameters for the pulse period and period derivative. We fit models in which the projected 
orbital radius (a x sin i) was a separate parameter for each of the 3 observations and also a model 
in which a x sin i was represented by a single parameter. 

In Table 3, we give the results of the fit to a model which includes separate orbital phase, 
pulse period, and pulse period derivative parameters for each of the three observations, but 
only global parameters for the projected orbital radius and for the orbital eccentricity. The 
pulse period and orbital phase zero are referenced to the solar system barycenter. We obtain a 
value for the eccentricity that is not statistically significant. The results presented in Table 3 
and Figure 1 are based on data from the energy intervals 7-14 keV, 6-24 keV, and 2-9 keV for 
the 1987, 1988, and 1989 observations, respectively. 

The value of 53.4875 ± 0.0009 it- sec that we obtain for a x sin i is consistent with the 
value of 53.46 ± 0.03 it-sec found by Primini, Rappaport, and Joss (1977). 

The residuals in the fits to the pulse arrival time delays are dominated by systematic 
effects. We arrive at this conclusion, in part, by noting that the residuals are very similar to 
those obtained when the analysis is performed on data from other energy ranges. Possible 
reasons for such systematics are (i) our lack of a complete model for the Ginga clock or (ii) 
effects that are actually due to some property of SMC X-l such as the spin behavior of the 
neutron star, small amplitude scattering of the X-ray pulses, changes in the intrinsic pulse 
profile, and so forth. Most of the evidence, however, points to the Ginga clock as being the 
source of the few-millisecond systematic residuals. 

Our determination of orbital phase zero is shown together with previous orbital phase 
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determinations for SMC X-l in Figure 2. The solid curve is the best-fit to a quadratic function. 
The coefficient of the quadratic term of the fitting function yields P orb / Port, = (-3.34 ± 0.02) x 
10-6 yr*l. This represents highly significant evidence for orbital decay. This value is about a 
factor of two larger than the observed orbital decay rate for Cen X-3. 

The orbital and other system parameters that we have determined for the SMC X-l system 
are summarized in Table 4. The values for a x sin i and P ort/^orb discussed above are included 
in the Table, as is the limit on orbital eccentricity of e < 7 x 10' 5 (2a confidence limit). This is 
presently the lowest limit on eccentricity for any of the known X-ray pulsars (true even for Her 
X-l?). The coefficients for the quadratic function that gives the time of the Nth eclipse in SMC 
X-l are also listed in Table 4. The value of Ppuise/^pulse = ' 5-37 x 10* 4 yr" 1 cited in Table 4 
is a fit to the long-term secular decrease in pulse period over the past 1 8 years. Note that the 
measured values of Ppuise/Ppulse determined for the three Ginga observations (Table 3) are 
each in reasonable accord with the long-term spin up trend. This is in sharp contrast to the 
situation in massive X-ray binaries where the accretion takes place via a stellar wind (e.g., 
4U0900-40; references). 

Values for the neutron star mass, companion star mass and radius, and orbital inclination, 
along with the corresponding uncertainties are also given in Table 4. These are derived from a 
Monte Carlo analysis which incorporates the measured system parameters, including the 
eclipse duration (see, e.g.. Joss & Rappaport 1984; Nagase 1989, and §4.1 for more details 
and references). 

3.2 Pulse Profiles 

We used the best-fit orbital parameters to produce pulse profiles for selected time intervals 
from each of the three observations. Pulse profiles for 6 different energy intervals are shown 
in Figure 3. All of the 62 ms resolution data from the 4-day 1989 July-August observation 
have been used to construct these profiles. In contrast to the pulse profiles from a number of 
other X-ray pulsars (see, e.g., Levine et al. 1987, Nagase 1989, Levine et^al. 1991) the shape 
of the SMC X-l profile is nearly independent of X-ray energy over the range 2-37 keV. 

Sample pulse profiles from the 1987, 1988, and 1989 observations are shown in Figure 4. 
While the pulse profile remains qualitatively unchanged (e.g., smoothly varying and double 
peaked), the detailed shape definitely changes with time time, especially the relative amplitude 
of the smaller of the two peaks. In particular, note that the smaller peak varied significantly 
over an interval of only 2 days (?) during the 1988 observation. On the other hand, it is 
apparent from the small residuals in our fits to the pulse arrival times (rms scatter of -2 ms) that 
the phases of the two peaks in the pulse profile remain highly stable. 

Finally, in Figure 5 we show a pulse profile for SMC X-l that was constructed from 7.8 


7 


ms resolution data (during the 1989 observation). The only additional structure that is revealed 
which is not seen in the profiles constructed with 62 ms resolution data (Figure 3) is the small, 
but statistically significant, feature on the rising protion of the smaller peak at pulse phase 
-0.95. Thus, the generally smooth nature of the SMC X-l pulse profile is in contrast to the 
profiles of some other X-ray pulsars, especially those with longer pulse periods (e.g., 
4U0900-40 and LMC X-4; McClintock et al. 1976, Levine et al. 1991). 

Paragraph about the shape of the pulse profiles: (i) not sinusoidal even thought they are 
"smooth", (ii) the two peaks are not 1 80° apart (so no simple centered dipole can explain the 
pulse profiles). 

4. Discussion 

4.1 Orbital Period Changes 

Observed changes in the orbital periods of massive X-ray binaries, including SMC X-l, 
are summarized in Table 5. Possible causes for changes in the orbital period of a binary 
system include mass transfer or mass loss, the gravitational influence of a third body, 
gravitational radiation, and tidal interactions. Kelley et al. (1983b) showed that conservative 
mass transfer was too small to explain the observed orbital period derivative of Cen X-3, and 
that mass loss from the binary system could explain the magnitude of the period derivative only 
if implausible assumptions about the rate of the mass loss and its specific angular momentum 
were adopted. Gravitational radiation is certainly negligible in these systems ( Port/P orb ~ *1 x 
10*11 yri). Furthermore, there is no evidence for a nearby third body. Tidal torques, on the 
other hand, could plausibly produce the observed negative period derivatives in both SMC X-l 
and Cen X-3; however, this requires that the companions of these systems are currently 
rotating more slowly than their orbits. 

Previously, the orbital decay of Cen X-3 and the limits on orbital period changes in LMC 
X-4 have been interpreted (Kelley et al. 1983; Levine et al. 1991) in the framework of the weak 
friction model (Darwin 1879, Kopal 1959, Alexander 1973). We briefly review that analysis 
and apply it to the observed decay of the SMC X-l orbit. In our analysis we make the 
simplifying assumption that the effects of mass transfer or loss can be neglected, so that the 
total angular momentum is constant and so that tidal interaction is the only mechanism for 
exchanging angular momentum between the rotation of the companion and the orbit. In this 
case, the total mechanical energy (orbital and rotational) decreases due to tidal friction unless 
the rotation frequency of the companion is synchronized to the orbit. 

We define a dimensionless orbital separation s = a/ty. total system angular momentum 
(orbital and rotational) j = J/(Ico 0 ) > and total system energy £ = E/( 1/2 Icdq 2 ), where r 0 is the 
orbital separation such that the orbit would have a moment of inertia equal to that of the 
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companion (r 0 2 = I/n), coq is the corresponding Keplerian angular velocity ((Dq 2 = GM T /r 0 3 ), I 
is the moment of inertia of the companion, |i. = M x M<./Mt is the reduced mass of the binary, 
Mj is the total system mass, and G is the gravitation constant (see also Darwin 1879, 

Counselman 1973). 

The total mechanical energy is then given by the following expressions: 

e = (j - s 1/2 ) 2 - s'* 0 a ) 

e = £2 2 s -3 . s-i , (lb) 

where £2 = cDc/cok, and (0 C and 0)^ are the rotation frequency of the companion and the orbital 
frequency, respectively. As discussed above, we have assumed that the total angular 
momentum remains constant as the binary system evolves tidally. Equation (la) can be 
evaluated to obtain the total energy as a function of s for different fixed values of the total 
angular momentum. The results are plotted in Figure 6. The value of £2 is also important in 
classifying the tidal evolution. Equation (lb) is plotted in Figure 6 for different values of £2 to 
illustrate the correspondence of £2 with other system parameters. Each point in the (s, £) plane 
(Figure 6) corresponds, in general, to two values of £2, and, hence, to two values of j . We 
assume, however, that the rotation of the companion is prograde (i.e., that £2 > 0); this implies 
that s < j 2 . Therefore, the curves of constant; in Figure 6 are terminated at s = j 2 . 

In this simplified analysis, the tidal evolution may be stable or unstable (Darwin 1879, 
Counselman 1973, and Hut 1981). In the former case, the end point of the evolution is an 
equilibrium configuration characterized by corotation of the companion with the orbit. In the 
latter case there is no equilibrium configuration, the system cannot evolve to achieve corotation, 
and the system will eventually undergo catastrophic- orbital decay. The dimensionless total 
angular momentum and orbital separation determine whether the evolution is stable. The 
classification of the courses of tidal evolution according to these quantities is summarized in 
Table 6. 

We have used a Monte Carlo error propagation technique (see e.g.. Joss and Rappaport 
1984, Nagase 1989) to indicate regions of the (s, e) plane in which the SMC X-l system is 
likely to lie. Briefly, the ingredients for such an analysis are the projected semimajor axes of 
the orbits of the neutron star and companion, the X-ray eclipse angle, the moment of inertia and 
rotation rate of the companion star, and a parameter (3 which indicates the degree to which the 
companion star fills its Roche lobe. Since we have no direct accurate observations which yield 
co c , we simply assume £2 to be uniformly distributed over the interval 0-1. In the context of 
the tidal interaction scenario, the decay of the orbit indicates that £2 < 1 . The assumptions and 
other input to the Monte Carlo system parameter analysis are discussed in Joss and Rappaport 
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(1984) and Nagase (1989). The moments of inertia of massive stars (both main-sequence and 
evolved stars) are based on unpublished stellar evolution calculations by Podsiadlowski 
(1990). 

The analysis of the SMC X-l system yields a range of solutions comprising cases in 
which both stable and unstable orbital decay are allowed (see Fig. 6). The region of the (s, e) 
plane occupied by SMC X-l in which stable orbital decay is allowed is, however, quite small. 
This region is characterized by 1.83 < s < 2.19 and 0.95 < £2 < 0.97. i.e., quite close to 
synchronization. All regions with Q < 0.95 correspond to systems with unstable orbital decay. 
We note that the region of the plane occupied by the SMC X-l binary system is very similar to 
that for LMC X-4 (Levine et al. 1991), except that in the former case the region with Cl > 1 has 
been eliminated. 

In the context of the tidal evolution scenario discussed above, binaries containing neutron 
stars are expected to be quite eccentric just after the supernova which created the neutron star 
(see, e.g., Lecar, Wheeler, and McKee 1976; van den Heuvel 1983). SMC X-l, Cen X-3, 
and LMC X-4 all have small eccentricities, which indicates that the orbits of these systems have 
been circularized. Hut (1981) has shown that, for sufficiently small orbital separations, the 
circularization time scale is shorter than the synchronization time scale (cf. Lecar, Wheeler, and 
McKee 1976). Thus, current theory allows the possibility that systems which have circularized 
by tidal interactions may not yet have synchronized. The observation of a decreasing orbital 
period in SMC X-l and Cen X-3 may therefore be the result of tidal torques from a companion 
star that is rotating less rapidly than the orbit because there has been insufficient time to achieve 
synchronization since the supernova which created the neutron star. The synchronization time 
scale becomes indefinitely long as s decreases toward V3. Indeed, if s < V3 or j < j CT then the 
system evolves unstably away from synchronization as- it decays. 

The above analysis, however, neglects the effects of the expansion of the companion due 
to its nuclear evolution, and the concomitant secular increase in mass transfer rate. These 
effects imply that catastrophic orbital decay is a necessary outcome, regardless of the stability 
of the tidal evolution predicted in the above analysis. As the companion star evolves, it will 
ultimately fill its Roche lobe (see §4.2). In turn, as the Roche lobe descends into the dense 
atmosphere of the companion, unstable mass transfer (on a thermal timescale of ~10^ yr) will 
commence. Since the neutron star cannot accrete matter at rates above the Eddington limit 
(~10'8 Mq yr 1 ) most of the matter transferred through the inner Lagrange point must be 
ejected from the binary system. If the specific angular momentum of the ejected matter is 
approximately that of the neutron star, then the orbital separation will shrink rapidly. Such a 
situation probably leads to a phase where the neutron star spirals into the envelope of the 
companion possibly leading to the formation of a Thome-Zytkow object (Thome & Zytkow 
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197X; Cannon & Eggleton 1992). However, while the SMC X-l system may be approaching 

this final phase of mass transfer, it is clearly not there yet. 

The expansion of the companion star due to its nuclear evolution can also drive tidal 
decay. As the companion evolves, it expands and its moment of inertia increases. The 
consequent decrease in its rotation rate drives the star away from synchronous rotation. The 
resulting tidal torque transfers angular momentum from the orbit to the rotation of the 
companion, thus leading to orbital decay. Below we estimate the evolutionary expansion rate 
of the companion that is necessary to explain the orbital decay of SMC X-l, and we comment 
on the plausibility of such an expansion rate. 

In the context of the weak friction model, the tidal torques on the companion star, N c , and 


on the orbit, N orb , may be expressed as 


N = -r = 

c dt 


daw) 


dU _ , (0> c - 1 (2b) 

orb dt dt X 

where x parameterizes the time scale of the dissipative tidal interaction. After Kepler s third 
law is initially used to eliminate the dependence on the orbital separation, these two equations 
yield the following differential equation that describes the departure from synchronous rotation 
as the companion star expands: 

d(ov -w) dln(I) (1 ~ 3I/M-a 2 )(<» K -® c > . (3) 

dt c dt x 

The first term on the right hand side describes the effects of the expansion of the companion, 
while the second term describes the effect of the transfer of angular momentum between the 
rotation of the companion and the orbit. This equation restates, that in the absence of 
expansion of the companion, and for sufficiently large separations (i.e., s 2 /3 = ^ia 2 /3I > 1), the 
system tends to evolve toward a synchronous state. When s 2 /3 = |ia 2 /3I approaches unity, the 
second term on the right in equation (3) decreases in magnitude. The first term may, at some 
point, become larger in magnitude than the second term and the expansion of the companion 
will then drive orbital decay. If JJ.a 2 /3I becomes less than unity, the expansion of the 
companion and the instability of the tidal evolution can both drive the orbital decay. 

Equation (2b) may be rewritten with the help of Kepler's third law to yield the following 

expression for the orbital period derivative: 
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31 (“c'V 


(4) 


P 2 T 

r orb [lad^ x 

The above equations are sufficient to estimate the expansion rate of the companion which 
is required to explain the observed orbital period derivative. To accomplish this, we assume 
that jxa 2 /3I is sufficiently close to unity that the second term on the right side of equation (3) 
can be neglected in comparison with the first term. When this condition applies for a length of 
time At, we can estimate the departure from synchronism as 


®c - "k = 


(Kw-Ok) 

dt 


At 


-(0 At 

C 


dln(I) 

dt 


(5) 


Substitution of this expression into equation (4) immediately yields 

^ ^Atdlnffl (6) 

P orb x ^ 

The measured value P orb /P or b ~ 3 x 10*6 yr' 1 may then be seen to be consistent with the 
following choice of parameter values: co c /co K “ 1, dln(I)/dt = 5x 10" 7 yr 1 (Podsialowski 
1990), and At/t ~ 6. This estimate is only valid if p.a 2 /3I is sufficiently close to unity since it 
assumes that 10,-0*! continues to grow for a time At - 6x. This will be true only if ll-3IAia 2 l 
< 1/6 or, equivalently, 1.6 < s < 1.9. Otherwise, the second term in equation (3) grows to 
sufficient magnitude to limit the increase of Io c -g) K I. One can also see that faster expansion 
rates of the companion would be able to explain the observed orbital decay rate without 
requiring s 2 /3 = }ia 2 /3I to be as close to unity. By comparison, the Monte Carlo analysis 

described above indicates that 1.6 < s < 2.2 for SMC X-l. 

Core hydrogen-burning stars with masses greater than -15 M© have dln(I)/dt > 5 x 10' 7 
yr -1 (Podsialowski 1990). However, it has been noted for some time (Hutchings et al. 1979, 
Rappaport & Joss 1983) that the masses, radii, and luminosities of the companions in massive 
X-ray binaries are not quite consistent with the properties computed for isolated stars with no 
strong stellar winds (Iben 198x, Podsialowski 1990). Specifically, a number of the massive 
companion stars appear to be overluminous and have larger radii for their mass (possibly as a 
result of substantial mass loss via a stellar wind; see, e.g., Conti 1978, Hutchings et al. 1979, 
Rappaport & Joss 1983). It is difficult, therefore, to determine the present expansion rate of 
the companion star in SMC X-l by comparison to the model calculations. As a consequence, 
we can say only that dln(I)/dt > 5 x 10* 7 yr' 1 is a plausible rate of evolution for the moment of 

inertia of the companion star, though it is hardly a secure estimate. 

We conclude that expansion of the companion cannot be ruled out as the driver of orbital 
decay, i.e., this scenario is plausibly consistent with our knowledge of the SMC X-l binary 
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system parameters and the evolutionary state of the companion star, but is otherwise unproven. 

4.2 Active Lifetime of Massive X-Rav Binaries 

SMC X-l and Cen X-3 are the two massive X-ray binaries for which both the orbital 
decay rate and the pulsar spinup time scale are known. The orbital decay time scales for the 
two systems are 3x10^ years and 5x10^ years, respectively; the corresponding spinup ti me* 
scales are 2000 and 5000 years. In each case, the ratio of these time scales is -100. Thus^the 
orbital decay time scales are about 10 times shorter than the nuclear evolution time of the 
massive companion stars in these systems, ImSS? are about two orders of magnitude longer than 
the spinup time scales. It is instructive to try to reconcile these apparendy disparate time scales. 

To address these questions, we have constructed a very simply model to explore the 
evolution of the spin of the neutron stars in these systems. Consider a system consisting of a 
neutron star in orbit about a massive stellar companion. For a description of the earlier phases 
of the evolutionary history of such a system see, e.g., van den Heuvel (1983, 1986, 1990). 
We begin our look at the system when the neutron star is orbiting its unevolved massive 
companion which substantially underfills its Roche lobe. In the -few x 10 6 years until the 
companion evolves to fill its Roche lobe, the neutron star will have spun down from its 
presumably rapid rotation rate at birth to a period of the order of a second (see, e.g., Lyne and 
Manchester review article). Interactions of the rotating pulsar with the stellar wind of the 
companion will spin the pulsar down even further (references on the propeller spin-down 
mechanism). At some point the neutron star will begin to accrete matter from the stellar wind at 
a relatively slow rate (e.g., — 10"^ 2 Mg yr*^). After a relatively brief interval of such accretion 
the rotation rate of the neutron star will attain an equilibrium spin period given by 



M 

- 3/7 

B 

P - 0.9 
«q 

10’ 8 M Q yr'' 


10 12 Gauss 


(Lamb, Pethick, & Pines 1972, Rappaport & Joss 1977), which turns out to be -18-150 s for 
a surface B field of 10 12 Gauss and values of M in the range of -10' 1 1 to 10' 13 Mq yr l . 

Eventually, the companion star expands as it evolves and begins to fill a more substantial 
fraction of its Roche lobe. As the companion fills a large fraction of its Roche lobe, the tidal 
torques discussed in §4.1 become effective and the orbit will start to decay, provided that the 
companion is rotating more slowly than the orbit. We assume that from this point on, the 
orbital decay time scale is shorter than the evolutionary time scale of the companion (but see 
§4.1). Within a small number of orbital decay time scales, the Roche lobe begins to descend 
into the dense atmosphere of the companion and the final stages of the luminous X-ray phase 
are under way. 
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Following the treatment of Verbunt and Rappaport (1988), we model the orbital evolution 
by the following equation constructed from a simple dimensional argument: 

- r/H 

( 8 ) 


dr 

dt 


-Re 


R s e 


ther 


where r = R^-Rg. Rl * s t * ie mean radius of the Roche lobe, and Rs is an arbitrary reference 
radius near the photosphere of the companion where the gas density is sufficiently high to 
sustain very large mass transfer rates. The other constants in eq. (6) are x ev , the time scale 
over which stable evolution of the binary orbit proceeded (e.g., via tidal interaction) before the 
onset of the unstable mass transfer, x^ is the time scale on which the thermal-timescale mass 
transfer proceeds as r — * 0, and H is the atmospheric scale height of the companion. Briefly, 
the second term on the right hand side of eq. (8) is derived from the approximation that Rl/R-l 
is proportional to M/M (see Verbunt and Rappaport 1988) and the assumption that the mass 
transfer rate grows exponentially with the descent of the Roche lobe into the atmosphere. The 
minus sign associated with this term follows from the fact that the orbit shrinks upon 
conservative mass transfer from the more massive to the less massive star, i.e., the transfer 
will not be stable. We note that even during the final stages of the unstable transfer, when M 
greatly exceeds the Eddington limit, the matter that is necessarily ejected from the system 
should carry away a specific angular momentum approximately equal to that of the neutron 
star, in which case the sign of the proportionality in equation (8) should remain unchanged. 

As shown in Verbunt and Rappaport (1988) the characteristic time scale for the growth of 
this runaway process is x — Hx ev /R$, i.e., the time for the tidally driven orbital decay to move 
the Roche lobe through an exponential scale height. For example, in the case of SMC X-l, if 
H is taken to be the pressure scale height at the photosphere (i.e., ~ 0.03 Rq) then x — 1000 
years, a number that is comparable to the pulse period spinup time. 

We can quantify the above discussion by integration of eq. (8) starting with the Roche 
lobe relatively far from the dense atmosphere of the companion, e.g., r > 100 H. As an 
illustr ati ve example, we take x ev = 3 x 10^ yr (see Table 4) and x^j. to be 10 4 years (the results 
are insensitive to this latter choice). At each step in the integration we compute the mass 
transfer rate from M/M “ exp(-r/H)/ x^r The mass transfer rate can, in turn, be used to 
follow the spin evolution of the neutron star. The equation that we integrated to study the spin 
evolution of the neutron star is 


dP 


spin ^ 

dt 


1.8 x 10' 4 P 


spin 


M 

6/7 

B 

1 O' 8 M 0 yr 1 


10 12 Gauss 

- 


in 


-1 


s yr 


( 9 ) 


(Rappaport & Joss 1976, Ghosh & Lamb 1979). As an initial condition, we assume that when 
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the Roche lobe of the companion is far from its atmosphere, there is still a small steady 

accretion rate from the stellar wind (e.g., 10' 12 M© y r ')- 

The results of such an integration are shown in Figure 7 for an assumed atmospheric scale 

height of 0.05 R©. The evolution with time of the quantity r = Rf R s ( in umts of is 
shown in Fig. 7a and the corresponding mass loss rate from the companion is given in Fig. 7b. 
The pulse period history along with the instantaneous value of the equilibrium pulse period are 
given in Figure 7c. Finally, the spinup time scale for the neutron star as a function of time is 
shown in Fig. 7d. Note that there is a specific time (i.e., at t = -16,000 yr) before the neutron 
star is engulfed in the atmosphere of the giant, when the mass transfer rate is M - 10 8 M© 
yr 1 , the spin period is -0.7 s, and the spinup time scale is P/P - 3000 years, all in agreement 
with the observed quantities. We find that for either substantially larger or smaller values of 
the atmospheric scale height, the above three parameters do not match the observed quantities 
simultaneously. The scale height that yields the 'best fit’ is 0.05 R© which is a factor of -3 
larger than that calculated for the atmosphere of a massive star of the type found in SMC X- 1 
(see, e.g., Kurucz 1977). We note, however, that there are significant uncertainties in the 
equation governing the spin up of the neutron star. These include the properties of the neutron 
star, especially the value of its magnetic moment. We estimate that, given all of the 
uncertainties that we can identify, the corresponding range of allowed values of the scale height 
that could yield satisfactory agreement with the observed parameters at the present epoch 
includes 0.01 R© < H < 0.05 R©. These values are not too far from the standard value for a 
nominal stellar atmosphere and do not require the kind of ad hoc extended exponential 
atmosphere of the types invoked in X-ray eclipse studies of such sources as 4U0900-40 and 
4U 1538-52 (see, e.g., Lewis et al. 1992; Clark et al. 1986?). 

The authors acknowledge the Ginga team for their support during the observation of SMC 
X-l. Two of us (AL and SR) are especially grateful for the hospitality accorded us during our 
stay at ISAS. We thank Ph. Podsiadlowski for providing us with some of his stellar structure 
and evolution calculations prior to their publication. Thanks to Jonathan Woo and John 
Chandler. 
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Table 1 


Orbital Phase Measurements of SMC X-l 



N 

T^(MJED) 

Method 

Uhuru 

-481 

40963.99 ± 0.02 

eclipse ingress - egress timing 

optical 

-249 

41866.94 ±0.04 

optical light curve 

Copernicus 

-144 

42275.65 ± 0.04 

eclipse timing 

SAS-3 

0 

42836.1828 ± 0.0002 

pulse timing 

Ariel 5 

42 

42999.6567 ±0.0016 

pulse timing 

COS-B 

72 

43116.4448 ± 0.0022 

eclipse timing 

Ginga 

1055 

46942.47240 ± 0.00013 

pulse timing - present work 

Ginga 

1173 

47401.744469 ± 0.000011 

pulse timing - present work 

Ginga 

1260 

47740.35910 ± 0.00004 

pulse timing - present work 


N is the number of orbital cycles, counted from the eclipse referenced for the SAS-3 data. 

T,^ (MJED) is the orbital epoch, defined as the time when the mean longitude equals k / 2 . For a 
circular orbit this is the time of superior conjunction. Times of mid-eclipse derived from ingress and 
egress observations are also listed in this column. 
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Table 2 


Ginga Pulse Height Energy Channels 


Channel Number Rftngg BllS£ Height Energies 

(keV) 


0-2.3 
2.3 - 4.6 
4.6 - 6.9 
6.9 - 9.2 
9.2-11.5 
11.5-13.8 
13.8 - 16.1 
16.1 - 18.4 
18.4 - 23.0 
23.0 - 36.8 
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Table 3 


Orbital Fits 3 for the Three SMC X-l Observations 


Parameter 
a x sin i 

Yalus 

53.4875 (9) It-sec 

Eccentricity 

<0.00007 (2 a) 

T^ (eclipse center, barycentric) e 

1987 May 

1988 August 

1989 August 

MJED 46,942.47240 (13) 
MJED 47,401.744469 (11) 
MJED 47,740.35910 (4) 

Ppulse (barycentric) b @ Epoch 

1987 May 

1988 August 

1989 August 

0.710185091(32) s (@MJD 46941.72954) 

0.7101406731 (12) s (@MJD 47402.52220) 

0.709809864 (43) s (@MJD 47740.35608) 

Ppulse/^pulse (Gm§ a ) 

1987 May 

1988 August 

1989 August 

- 5.5 (7) x 10' 4 yr' 1 
-4.54 (2) x 10' 4 yr' 1 
-5.5 (2)x 10' 4 yr' 1 


a All uncertainties are la single-parameter confidence limits, except for the upper limit on orbital 
eccentricity which is 2a. The numbers in parentheses are the uncertainties in the last digits given. 

b? 

c? 

d? 

List the definitions of MJED and MJD in the caption 
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Table 4 


Parameters for the SMC X-l Binary System 


Parameter 

Value 

Uncertainty 

a x sin i 

53.4875 lt-s 

0.0009 

eccentricity 

< 0.00007 (2 o) 

0.000016 

Ppulse/Ppulse dong term) 

-5.37 x 10' 4 yr 1 

2 x 10' 6 

ao^ 

MJED 42,836.18277 

0.00020 

ai o» 

3.89229073 days 

5.6 x 10* 7 

a 2 ( b ) 

-6.92 x 10' 8 days 

0.04 x 10' 8 

Port/P orb 

-3.34 x 10' 6 yr 1 

0.02 x 10' 6 

Me® 

18 Mg 

2.0 

M x ( c ) 

1.1 Mq 

0.2 

^eclipse^ 

26.5-29° 

... 

inclination angled 

62° 

4° 

Rc (c) 

17 R 0 

1.5 


What is footnote a? 

b The time of the Nth eclipse can be written as: t N = ao + a t N + a 2 N^ 

c Based on a Monte Carlo error analysis (see, e.g.. Joss and Rappaport 1984, and Nagase 1989). 
d References in Joss and Rappaport 1984 
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Table 5 


Orbital Period Changes in Massive X-Ray Binaries 


Source 

PortAWb (y^) 

Reference 

Cen X-3 

Porb^orb = -1-* X 'O' 6 

Kelley et al. 1983b 

SMC X-l 

Porb^orb = -3-3x10-6 

Present Work 

LMC X-4 

IPorl/Porb 1 < 2 x 10-« 

Levine et al. 1991 

Her X-l 

Porb/Po,b“ -1-3x10-8 

Deeter et al. 1991 

4U1538-52 

^ort/^orb^ < 6 x 10 ^ 

Cominsky & Moraes 199 1 


22 


Table 6 


Constraints on Tidal Evolution 



&<lL 

> 1 

<3 

Unstable Decay 

Stable Growth 

^3, j <j C r 

Unstable Decay 

Not Possible 

^3, j > j cr 

Stable Decay 

Stable Growth 


The symbols are defined as follows: s is the dimensionless orbital separation, j is the 
dimensionless orbital angular momentum, and Q is the ratio of companion rotation frequency to 
orbital frequency. See text for definition of the unit of length and angular momentum. 
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Figure Captions 


Figure 1 : Doppler delay data for SMC X-l obtained from Ginga timing observations. The 
filled circles are the measured pulse arrival times with the effects of a constant pulse period 
removed. The solid curve is the best fit to a circular orbit (see text and Table 3). The residuals 
between the best fit and the data are shown as an "x" (note the expanded scale for the 

residuals). 

Figure 2 : Summary of mid-eclipse delay times for SMC X- 1 with respect to a constant orbital 
period of 3.89229073 days. The SAS-3, Ariel 5, and Ginga data points are taken from Kelley 
etal. (1983b), Ariel 5 reference, and the present work, respectively. The filled circles are the 
measured delays of mid-eclipse as determined from pulse timing analyses; the filled triangles 
are from X-ray eclipse timing measu re me nts ? w h ile the open square is from optical studies. The 
solid curve is the best fit of a quadratic function to the mid-eclipse delay times (see Table 4). 

Figure 3 : Pulse profiles of SMC X-l, in 6 energy bands, from the 1989 observations. 
Non-source background has been subtracted. The pulse profile is defined by 64 phase bins 
(-0.21-s time bins). For each pulse profile, a typical ± la error bar due to counting statistics is 
shown. Data obtained during the eclipses, large flares, or times of particularly large pulse 
profile variability were excluded from the construction of the pulse profiles. 

Figure 4: Sample pulse profiles from the 1987, 1988, and 1989 Ginga observations. Note 
that the relative amplitude of the smaller of the two peaks in the profile changes with time, 
including significant variations over an interval of only 2 days (?) during the 1988 observation. 

Figure 5: Pulse profiles for SMC X-l constructed from 7.8 ms-resolution data (during the 
1989 observation). Note the small, but statistically significant, feature on the rising protion of 
the smaller peak at pulse phase -0.95. The solid curves superposed on the histograms are the 
pulse profiles derived from the 62 ms-resolution data. 

Figure 6 : Binary system energy vs. orbital separation including the kinetic energy associated 
with the rotation of the companion star. The quantities E 0 = l/2Ico 0 2 and R 0 = ^V\i are scaling 
constants to cast the energy and orbital separation in dimensionless form (see text for details). 
The solid curves represent constant values of total angular momentum, while the dashed curves 
are for constant values of the parameter Q = co c /co K (see text). The curve of constant total 
angular momentum labeled "j cr " corresponds to the smallest value of j for which these curves 
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have a local minimum. The heavy contours represent 50% and 95% confidence limits on the 
measured SMC X-l system parameters, as determined from a Monte Carlo error propagation 
technique (see text for details). Because of the finite bin size in which the contours are 
calculated, they slightly overlap the non-physical region wherein E/E c < - s' 1 . 

Figure 7: Illustrative calculation of the evolutionary spin history of a neutron star in a massive 
binary system, (a) Difference between the Roche-lobe radius and the photosphenc radius in 
units of the atmospheric scale height (taken to be 0.06 Rq). (b) Rate of mass loss from the 
donor star, (c) Pulse period (solid curve) and equilibrium spin period (dashed curve), (d) 

S pinup timescale. 
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